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Abstract 

We investigate the width of the flux tube between heavy static quark charges. Using 
the gauge/gravity duality, we find the properties of the minimal connected surface related 
to the width of the bound state. We show that in the confining phase, the logarithmic 
broadening predicted by the effective string description and observed in lattice simula¬ 
tions is a generic property of all confining backgrounds. We also study the transverse 
fluctuations of the string connecting two static quarks in curved backgrounds. Our for¬ 
malism is applied to AdS space where we compute the expectation value of the square 
of transverse deviations of the string, a quantity related to the width. 


Contents 


1 Introduction 

2 Tube width from circular Wilson loops 

2.1 The setup. 

2.2 The connected minimal surface . 

2.3 Solutions for special backgrounds. 

2.4 Long flux tubes. 

3 Flux tube and connected minimal surfaces: A numerical treatment 

4 Fluctuations of the spatial string 

4.1 Static fluctuations. 

4.2 Time dependent fluctuations. 

4.2.1 The flat background . 

4.2.2 Curved backgrounds. 

4.3 The width of the fluctuation. 


Q 

B 

0 

0 

0 

0 


0 


| 12 | 

12 

11 

15 

H 

18 


5 Conclusions 

A An analytically solvable toy model 



21 


1 Introduction 

The confinement property of the strong force is of fundamental interest and while a 
lot of effort has been devoted to its study, our understanding of the non-perturbative 
dynamics involved is far from being complete. A way to study confinement is to use 
as probe a heavy bound state of a quark and an anti-quark. The fundamental quantity 
that one can compute using the pair is the expectation value of the rectangular Wilson 
loop. The shape of the loop is related to the fact that the probes are assumed to be 
very massive and that through time they stay at fixed positions. The two long legs of 
the Wilson loop correspond to the time direction, while its other two sides are along the 
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alignment of the pair representing the distance between the color charged particles. In a 
confining theory the Wilson loop gives the known area law, which is equivalent to a linear 
static potential in the inter-quark distance. The expectation value of the Wilson loop is 
therefore considered to be a criterion for confinement. At low temperatures, below the 
deconfinement transition temperature, string models predict a decrease of the potential 
and in particular the string tension decreases as temperature increases mmm, fact that 
has been also confirmed recently using the gauge/gravity duality pfj. 

There exists another quantity of almost equal importance that has not attracted 
the same attention as the potential of confining theories. It is the chromo-electric flux 
tube formed between the quark pair. Due to confinement, the chromo-electric field 
energy density is believed to be constrained in a tube-like region. In confining theories 
the tube has some interesting universal properties. Using a low-energy effective string 
description it has been shown that the flux tube width grows logarithmically with the size 
of the bound state, when measured at the middle of the pair [5]. This effect, sometimes 
called roughening, together with the appearance of the Luscher term signals the quantum 
string-like behavior of the tube joining the pair. The rough phase is described by strong 
fluctuations of the collective coordinates of the position of the string and the mean 
squared width of the flux tube diverges logarithmically when the inter-quark distance 
goes to infinity. The string like features of the confining flux tube have been studied in 
several gauge models, summarized in recent reviews, including wzm- The broadening 
is believed to occur due to the transverse fluctuations of gluonic strings that increase 
in magnitude as the quarks are separated. At zero temperature, the logarithmically 
divergent string width has been verified by several lattice computations with increasing 
accuracy mmm- These results have been extended by the study of the thickness of 
the tube in higher representations of the gauge groups. It has been found that the width 
of the flux tube of the bound state grows logarithmically due to the excitations of the 
underlying k-strings (Ej. Moreover, the width of the baryonic flux-tube junction from 
effective string theory has also been found to grow logarithmically ra. It is interesting 
to point out that logarithmic like behavior has been observed in several other similar 
models, in statistical physics and condensed matter physics, therefore hinting to universal 
properties. A more recent result is that the roughening undergoes a qualitative change at 
finite temperature slightly before the deconfinement phase transition, as the broadening 
apparently changes from the zero temperature logarithmic law to a linear law mmm 
as the phase transition is approached. All these properties of confinement can be studied 
holographically confirming the lattice theory results and with a hope that new unknown 
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properties will be revealed. 


There are more than one equivalent ways to measure the effective width of the flux 
tube. In lattice studies, the computation is done by considering two Wilson loops. One 
Wilson loop W(C) corresponds to the actual heavy QQ bound state that generates the 
flux tube, with one large time direction and one smaller spatial direction fixed to be the 
size of the bound state. The other, small probe Wilson loop P(c) is placed on a plane 
above the QQ Wilson loop plane in the center and measures the flux passing through it 
|S]. The chromoelectric field energy density £{x) is defined as 

(W(C)P(c)) - <iy(c))<p(c)) 

me)) 1 ' ' 

and recall that in holography it is specified by the minimal surface connecting the two 
Wilson loops. The mean square width of the flux tube is then defined as 

2 f dx±x 1 2 ± £(x) 
f dx±£(x) ’ 

where x± are the plane dimensions transverse to the direction of the heavy pair. 

Another way to calculate the width of the flux tube is through the time-dependent 
transverse fluctuations of the string connecting the static quarks with appropriate Dirich- 
let boundary conditions for the spatial direction. In that case the expectation value of 
the square transverse deviations of the string has been studied extensively in mm and 
is given by 



w 2 {x) = (x(x,t) 2 ) =lim(x(x,t)x(x'= x + e,t'= t + e)) 

= lim G(x,t]x'= x + e,t'= t + e)) (1.3) 

for each independent, transverse, fluctuating degree of freedom y, where G(x, t; x', t') is 
the appropriate to the geometry Green function. For fundamental strings and in flat 
space, the width of the flux tube increases logarithmically with the separation distance 
of the charges. 

Here we extend these studies in the context of gauge/gravity dualities [16, ITT] . An 
effort in this direction has been already made in the case of strings in the fundamental 
representation with a confining AdS soliton background m and with k-strings in hard 
wall AdS G3]Q In our work we approach the problem from two different but equivalent 

1 Other related constructions of connected minimal surfaces aiming to study asymptotically the min¬ 

imal area appear in m- 
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ways. Finding the connected minimal surface between the probe Wilson loop and the 
Wilson loop corresponding to the heavy QQ pair, we use certain properties of the minimal 
surfaces in confining backgrounds and we prove that the logarithmic broadening is a 
property of any confining gauge/gravity duality. The approximations used for our generic 
proof are validated in the numerical section for the confining DA Witten model, where 
the exact solutions of the minimal surfaces are found. 

Then we look at the fluctuations of the gluonic string. Its static fluctuations can be 
worked out in a generic framework. The time dependent fluctuations are more interesting 
and challenging, where the generic formulas are again obtained. We begin by reviewing 
the time dependent fluctuations in flat space where we reproduce the logarithmic broad¬ 
ening. Then we move on to curved backgrounds and focus on the simplest conformal 
AdS case, where even though the meaning of flux tube is obscure (since the gauge theory 
is non-confining), an analogous quantity to the width w 2 can be defined. We describe the 
methodology we follow and we find numerical evidence that the quantity we call ’width’ 
in this space depends at least on one logarithmical term and is decreasing for small bound 
states. A toy model where the width of a string in the curved space may be obtained 
analytically is described in the appendix. 

Our paper is organized as follows. In section 2, we compute the width of the QQ flux 
tube in a generic confining gauge/gravity duality, using the connected minimal surface 
between two Wilson loops, one of larger and another of smaller radius. In section 3, we 
construct the exact numerical solution of the connected minimal surface in the D4-Witten 
background. In section 4, we consider fluctuations of a spatial string in the curved spaces, 
where we particularly focus on the AdS space. Finally, in section 5 we summarize and 
discuss our results. 


2 Tube width from circular Wilson loops 

In this section we compute the width of the tube using probe Wilson loops. To simplify 
the setup we consider instead of rectangular, circular Wilson loops. This approximation 
is valid for the study of the flux broadening since in the limits where the radius of the 
large loop is much larger than the size of the small loop and the separation distance, the 
shape of the boundary loops is not expected to play the leading role in the calculation 
of the width w 2 . 
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2.1 The setup 


In our analysis we use the generic metric 

ds 2 = -g 00 (u)dxl + gu(u)dx . 2 + g uu (u)du 2 + dX 5 , (2.1) 

where we assume that it is diagonal and isotropic, so all the space metric elements are 
equal. The X 5 is the internal space where all our configurations are kept localized. As 
we consider circular Wilson loops, it is advisable to make a change of coordinates on the 
space-like components of the metric, x\ = r cos 0 , x 2 = r sin 0. giving 

ds — — goodx 0 T gw (dr T v dO T dx^j T 9uud r u> T dX^ . (2.2) 

The boundary of the metric is taken to be at u = 0, where lim u _>o 9n( u ) = oo. 

Large part of our analysis aims to produce the universal behavior by using certain 
properties of the minimal surfaces in the confining metric. At the end we focus on a 
particular confining theory working numerically. 

2.2 The connected minimal surface 

We consider two circular Wilson loops at the boundary of the space with radii Ri and R 2 
(R\ < R 2 ), with the same center, but separated by a distance L along the x 3 direction. 
The string world-sheet is parametrized by 


6 = r , x 3 = a , r = r(a) , u = u(a) , 


(2.3) 


where the static gauge is chosen along 6 and x 3 . 
The Nambu-Goto action becomes 


S = — [ dadrgnr^ /l + r' 2 + := — J dadrgnr\fD 


2n 


9n 


2n 


(2.4) 


with T the string tension. Time does not appear explicitly in the action so the Hamilto¬ 
nian is a constant of motion, say —c, specified by 


c = 


9nr 

Vd' 


(2.5) 


that translates into the constraint equation 


! + r' 2 + ^V 2 - ^ = 0. 


9 li 


c 


( 2 . 6 ) 
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Regarding the dynamical coordinates r and u, using eq. (12.51) and after some algebra we 
arrive at the relatively simple equations 

,2 


r"-^ = 0 


// 9uu , / 2 1 Q ( Qvm \ r 2 g u d u gu 


w- - \-u -d u 

9n 2 \gn 


= 0 . 


(2.7) 

( 2 . 8 ) 


Hence the system of equations that needs to be solved consists of eqs. (12.61) . (12.71) and 
(12.81) . For particular class of backgrounds it can be solved analytically. 


2.3 Solutions for special backgrounds 


Let us consider the special case g uu = gu which leads to simplifications. Equations (12.81) 
and (12.61) then reduce to 


// r 2 g u d u g u 
u - - -= 0 , 


l + r /2 + u /2 -^ = 0. 


Multiplying (12.71) by r and (12.91) by u and adding it to (12.101) we obtain 

i( r y + i(«y - —( 4 ^ + M a„( 5 ? 1 )) + 1 = 0. 


(2.9) 

( 2 . 10 ) 


( 2 . 11 ) 


This equation can be integrated analytically only when gu = c/u 2 , corresponding to 
the AdS case, already done so in EH- Other simplifications for the choice of the metric 
elements which lead the differential equations to analytic solutions do not seem possible. 
Therefore, in order to obtain analytic solutions for generic backgrounds one needs to go 
to certain limits in our configuration. Thus, in the next section we consider long tubes 
in confining theories. 


2.4 Long flux tubes 

Concerning the limit of long flux tubes in confining backgrounds, one finds that sev¬ 
eral simplifications are possible. These come from the properties of minimal surfaces in 
confining spaces which are described in detail here. 

The string profile along the radius r of the tube lies very close to the bulk scale Uk 
defining the tip of the confining geometry. The projection of the surface to (r, u) has 
almost a rectangular 77 shape departing from the boundary and entering into the bulk 


6 













steeply. Departing from the boundary and the small Wilson loop of radius Ri, the surface 
enters steeply in the bulk attempting to minimize its area and to shrink the radius of the 
tube, thus moving on the left of the r axis. Then being at the bottom of the geometry 
the radius of the surface increases as we move to the direction of the larger Wilson loop. 
It extends for a distance which is almost equal to the separation of the radii on the 
boundary R 2 — Ri, and then just before the radius reaches the i? 2 value returns steeply 
to the boundary. 

The projection of the world-sheet along the (u, x 3 ) of the two Wilson loops is expected 
to lead again to an inverted II shape with a mostly rounded corner on the side of the 
small loop. The minimal surface at the side of the large radius Wilson loop goes deeper in 
the bulk to reach the smaller area elements, than the side that is closer to the small loop. 
The turning point of the string along the direction x 3 is not in the middle of the string but 
lies closer to the side of the large Wilson loop. A large portion of the string world-sheet 
lies around the turning point and it is this region that gives the dominant contribution 
to the minimal surface. The qualitative picture of the minimal surface solution described 
here is the same for any background and will be later verified numerically, by finding the 
particular solution of the surface in the D4 -Witten model. 

To proceed with the calculation let us define the two quantities 

h(u) := 4. /(«) ~ — , (2.12) 

011 

that allow us to write the system of equations in the more compact form 

r" — hr = 0 , 

2u" + {u') 2 d u (In/) — r 2 ^j^ = 0 , (2.13) 

(r') 2 + f(u') 2 = hr 2 — 1. 

At the bottom of the world-sheet where u = uq and u' — 0 we can determine its shape 
by solving the system 


r" — hr = 0 , ( r'y + 1 = hr 2 , 


(2.14) 


yielding to 


r = 


gSvT(ci+<r) ^g—2 svT(<t+ci) _|_ 


2 h 


(2.15) 


where C\ can be taken to zero without loss of generality and we choose the sign s = — 1. 
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The radius of the flux tube at its endpoints is specified by the boundary conditions 


r(0) = R u r(L) = R 2 . (2.16) 

Since the string world-sheet dives directly into the bulk and extends down to the tip of 
the geometry from where it takes its main contributions, it is natural to assume that 
the radius of the flux tube in the bulk around a = 0 and a = L approaches the radius 
of the string on the boundary And since the string lies on the tip of the horizon we 
can apply these boundary conditions to the function r{a) obtained in (12.151) . to get the 
leading contributions 

*-H i+ 9’ R ’=K e ^ L+ y w )' <2i7 > 

where h now is computed for u ~ i// c , close to the tip of the confining geometry. Notice 
that for large L a logarithmic dependence of the length on the ratio of the radii appears 
as 

(2 , 18) 

Minimizing the action and using the Virasoro constraint we obtain 

S=-V\^l f L dar \ (2,19) 

c Jo 

Integrating the action and using (12.171) and (12.181) we obtain 

s = ^ (2= + Rl - Rl + i(l - r 2 ^)) + % - «?) , (2-20) 

where we have used that u 0 ~ i// c and taken R 2 A> R\ . Notice that for the string 
tension we have substituted gu(uk)/ 2. The energy of the string obtained for confining 
backgrounds, in this limit, approaches the one obtained when the string is placed in flat 
space. 

Inserting the area of the minimal connected surface in (11.11) for the transverse direc¬ 
tions, we obtain 

w 2 ~ log ^ . (2-21) 

-Ri 

Evidently the width of the tube has a logarithmic dependence on the radius R 2 , i.e. on 
the inter-quark distance. Notice that the position of the small probe loop does not affect 
the result, as expected. 













3 Flux tube and connected minimal surfaces: A nu¬ 


merical treatment 


In this section we present a numerical treatment of the connected minimal surface between 
the Wilson loop associated with the heavy quark bound state and the probe Wilson loop. 
We remark that the tube-like surfaces that we plot do not correspond to the actual QCD 
flux-tubes, but result from the definition (11.11) . We obtain the exact minimal surface 
numerically, and plot the projections of the radial profile of the string u(a) and the cross 
section r(a) in terms of the distance a = x 3 . The findings of the numerical analysis show 
the validity of the analytical approximations made in the previous section. 


The theory we choose is the confining D4 Witten model m, with metric written in 
the form (12.211 


ds 2 = 


(— dt 2 + dr 2 + r 2 d0 2 + dx\ + f(u)dx \) + 
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( 7 / \ 3/4 7/3 

-) , f(u) = 1--J, R 3 = 7rg s N c l 3 s , (3.1) 

where dQ 2 is the metric of a unit S' 4 of radius R on which the string is localized. The 
spatial compactification is along x\ and the tip of the geometry is denoted by up-. The 
string tension in this model is l/(2na.') (n^/i?) 3 ^ 2 . 


The Wilson loops are placed in this background with the parametrization (12.311 with 
the boundary conditions of the previous section. Here we solve the system of equations 
(I2.13j) numerically, with no approximations. The input to the system is a choice for the 
quantities ro, uq, v!(gq ) = 0 and the value of r'(cro). The derivative of r at the turning 
point do is obtained by solving (12. 6 j) and using the values chosen already for the other 
parameters. The solution to the differential equations gives two functions of the form 
r(s) and u(s). 


The setup for the minimal surface is the following: The boundary of the space is at 
infinity as can be seen from (13.111 . We place the small Wilson loop with radius R\ at 
smaller distance from the origin along the direction of the tube x %, and the large Wilson 
loop with radius f? 2 at larger distance. The order of scales is R> > R\- The particular 
numerical values of the scales will be chosen to ease the presentation. The area of the 
minimal surfaces is regularized by either inserting a cut off close to the boundary or by 
substracting the infinite parts of the action mm- 

In Figured] the solution for u(r) is plotted. We observe the almost i7-like shape of 
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Figure 1: The profile of the minimal surface 
along the holographic direction u and its radius 
r. The shape of the curve resembles an inverted 
77, while the bottom of the string is moved to 
the left where the smaller loop is located. The 
minimum radius of the tube is inside the bulk 
and closer to the small loop. 


Figure 2: The minimal surface projected along 
u(x 3 ). For large values of X 3 where the larger 
Wilson loop is placed, the surface goes slightly 
deeper in the bulk, in order to minimize its area. 
The shape of the curve approaches again that of 
an inverted II, with the side closer to the small 
Wilson loop slightly rounded. 


r 



Figure 3: The radius of the string along the spatial extension of the tube. It is not a monotonic 
function and appears to have minimum close to the small Wilson loop. 

the string described in the previous section, where most of the world-sheet lies along the 
tip of the geometry. As the string enters into the bulk, the radius of the tube decreases 
at both ends. As expected, the smaller tube radius is observed on the left where the 
smaller Wilson loop is located. Notice that the minimum radius is inside the bulk and 
not on the boundary. 

In Figure[2]the function u(x 3 ) is shown. The rounded from one side, II shaped Wilson 
loop is presented. The minimal surface prefers to go deeper into the bulk to smaller area 
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elements at the side of large radii in order to minimize its area more effectively. That 
is why the minimum value along the holographic direction is closer to the large Wilson 
loop. Around this point the string lies almost flat for a long range. Moreover, this is the 
region from where the major contributions to the action come. 

In Figure |3] the solution for r(x 3 ) is presented, representing the radius of the tube 
along the X 3 direction. Starting from the small loop the radius decreases along the tube. 


The minimal value of the radius is closer to 
increases until it becomes equal to the one of 



0 


Figure 4: The two dimensional surface of the 
radius of the tube as it enters into the bulk. No¬ 
tice that the minimum of the radius of the tube 
occurs very close to the minimum holographic 
radial distance. The increase of the radius from 
its minimum value to the value of the radius R 2 
of the larger Wilson loop happens almost totally 
around the tip of the geometry u 


the small loop, after the turning point it 
the large tube. 



Figure 5: The minimal surface in the spatial 
dimensions. The small Wilson loop is placed 
at small X 3 and the large one on the other side. 
When the surface departs from the Wilson loops 
its radius is decreasing considerably. The spa¬ 
tial X 3 extension of the tube in the particular 
plot is five times enhanced for presentation pur¬ 
poses. 


In Figure [4] the two dimensional surface (r(u),£) of the tube entering the bulk is 
presented in X{ coordinates, used to bring the metric in the form (12.21) . The ^-dimension 
of the tube is not depicted. The inner funnel corresponds to the Wilson loop of radius 
R\ and as the surface enters into the bulk its radius becomes smaller. The inner funnel 
radius has reached already its smallest value before it reaches the tip of the geometry at 
the minimal holographic distance. The radius increases rapidly from its minimal value 
towards R 2 close to u 0 . Having approached a value close to R 2) the surface returns into 
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the bulk. The reason for this is that the main contribution to the area of the minimal 
surface comes from the region around uq. 

In Figure [5] the two dimensional surface (r(x 3 ). t) of the tube is shown in the Xi 
coordinates. Along the spatial direction x 3 the radius of the tube decreases until it reaches 
a minimal point and increases again when the loops are approached. The representative 
values we have chosen for the particular plots are c = 0.5, uo/uk = 1.005, tq — 1 that 
then fix i?i = 0.89, R 2 = 2.74 for the radii and Li(l) = 1.48, L 2 {2 ) = 2.17 for the 
position of the loops along the x 3 direction. 

4 Fluctuations of the spatial string 

In flat space the effective string description shows that the broadening of the flux tube 
at finite temperature, below the deconfinement phase transition, becomes linear. The 
linear broadening is supported also by Monte Carlo simulations in certain simple models 
m This means that there may be another order parameter characterizing confining 
gauge theories. From the point of view of string theory on the other hand it may be also 
used to distinguish confining backgrounds, as the expectation value of the Wilson loop 
does. The quantity w 2 in eq. (11.3j) is well defined not only for confining theories but for 
a larger class of backgrounds. 

The simplest object for which we can compute its fluctuations is a spatial string. Our 
bulk order parameter will therefore be 5x 2 (a 0 ), with x 2 one of the boundary, transverse 
degrees of freedom and cy, the world-sheet parameter specifying the middle of the string. 

4.1 Static fluctuations 

Let us consider the spatial string 


t = t , X\ = Xi(a) , x 2 = 0 , u = a 


(4.1) 


that leads to the action 



(4.2) 


and the equation of motion 



y/goopgiio-F 

9oo9n 


(4.3) 
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with g Uo := gn(a 0 ), g 00o := g 00 (o 0 ) and 


F : 


9oo9uu9u 


googn — 9ooo9iio 


(4.4) 


Perturbing 

X\ = Xi(a) + 8 x 1 (cr) , X 2 = 8 x 2 ( 0 -) , (4-5) 

one sees that the first order contributions cancel out [2fiJ while to second order the 
contribution of the fluctuations is 

c [ j (9oo9u U gu x 12 . 9 00 922 x , 2 \ (A 

S * = J da [ 2F3/2 + ~ 27¥ 5x2 ) (4 ' 6) 

and the resulting linearized equations of motion are 

d ( gSo9uu9u 5 _ Q ± ( 9oo922 5x ,\ =0 U7) 

do\ F 3 / 2 V ’ da l FI/ 2 V 1 j 

Since we care for the transverse fluctuations hr :2 and evidently the two types of fluctua¬ 
tions decouple from each other, we can set consistently the longitudinal part to zero. 

The two equations then that govern the string configuration and the fluctuations 


are 


r y/F 

aq(cr) - aq(<7 0 ) = y/go%gu 0 h, h(a\a 0 ) = / da - 

Ja 0 9 oo 9 n 


Sx 2 (a) - 8 x 2 (a 0 ) = -Ci/ 2 , 


r y/F 

h(cr-,a 0 )= / da - 

J (tq 900922 


(4.8) 


For spatially isotropic backgrounds gn = /y 22 and thus I\ — J 2 = / and .i'i(Vxo) is fixed 
by symmetry to be the middle point of the string. The solution is determined by the 
boundary conditions 


( = VSlloJoOci I(oo\<7 0 ) , Sx 2 (tr 0 ) = Ci/(oo;<7 0 ) , 


(4.9) 


that yield for our quantity of interest the 


PYnrpssinp 


8 x 2 ( 00 ) 


Cl 

2v / 5 , 00ofi , ii o 




(4.10) 


where r/ llo should be determined by substituting oq(L) from the first expression of (14.91) . 
The validity of the solution depends of course on the continuity conditions of the string 
in the turning point and should be checked for each background. 


Particular backgrounds 

The (14.81) can solved for two particular backgrounds, a conformal and a confining. 
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The metric elements for the AdS background are given by gu = g 2 2 = u 2 /R 2 and 
9uu = R 2 /u 2 , where R is the radius of the AdS space where we set it equal to the unit 
and we recover it only in the final result. Then 


/(cr;cr 0 ) = 


da- 


1 


'cro C T 2 y /( 0-4 - °o) 

The length L on the boundary is then given by 123 

2^r (f) 


L = 


*,r (J) 


y^r (|) 
(i) 


(4.11) 


(4.12) 


The metric elements for the confining ZM-brane background are —goo — gn — 922 = u 3 / 2 
and g uu = 1/ (u 3 / 2 /), where / = 1 — 'u|./u 3 , where we have set the radius R of the sphere 
equal to one and recover it for the final result. The length of the string is given by 


L 


2 a, 


. 3/2 


da 


'T 0 


V(r 3 o)(T 3 


—2 log (oo — it*;) + const. 


(4.13) 


where the constants containing terms from the boundary, are omitted. The turning point 
of the string is now determined by c 2 = v /goo 0 fi , iio and the analysis is similar to [2$, [21) , 
although in a different gauge. Then we obtain that 


(To — Uk (l + e 


(4.14) 


4.2 Time dependent fluctuations 


In this section we consider time-dependent fluctuating strings in curved backgrounds. 
The expectation value of the square transverse deviations of the string with periodic 
boundary conditions in Euclidean time, has been studied extensively in mm and is 
defined as explained in (11.31) 


w 2 (x) = (y(x, t) 2 ) = lim G(x , t; x' = x + e, t' = t + e)) , 
6 — 1-0 


(4.15) 


for each independent, transverse, fluctuating degree of freedom y, with G(x, t: x'. t') the 
appropriate to the geometry Green function. The coordinate x is some space-like world- 
sheet coordinate while t is the corresponding time-like coordinate. 


With time dependence built in as in [EE], the action of the perturbations is computed 
from (14.21) to give 



da 


9009 ^ 911^,2 , 900922 ^,2 

F 3/2 0X 1 + Jp ° X 2 


9n9uu 


5x\ 



(4.16) 
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leading to the equations of motion 


d 

da 

d 

da 


9oq9uu9ii 

_F 3 /2 


8 x[ 


900922 x ^, 
f 0 00 o 

Vf 2 


QnQuu 

~7T 


8 x i 


922^/F 

-0X2 

S'oo 


= 0, 


0 . 


(4.17) 

(4.18) 


We will restrict ourselves to the isotropic case gn = # 22 • The solutions for the string 
worldsheet are the ones we have already obtained while the solution for the perturbed 
string need to be reexamined due to the additional fluctuations. 


4.2.1 The flat background 


According to (11.31) we need to calculate the Green function of a Laplace equation on a two 
dimensional plane with one Dirichlet and one periodic boundary condition. The Green 
function of this cylinder with dimensions Lx f3 (/3 can be thought of as a Euclidean time) 
can be mapped to that of a torus with size L x f3 where L = 2L. The Green function of 
the Laplace equation with toroidal boundary conditions is given by 


G(iM) = WA + 4- 

2a/3L 27 ra 


E 

71—1 


COS 


(2nnx\ 1 
V L ) n(l- q n ) 


exp 


—2mrt 


L 


q n exp 


2 mrt 




(4.19) 

where q = e 2?7rT , r = i(3/L, K = /3/(12aL) + log 77 (r )/(77 a) and a is a string tension, not 
to be confused with the world-sheet order parameter. 


The cylinder Green function is then obtained from that of the torus by 


G c (x, t ; x ', t ') 


Gt(z — E ) — Gt(z z >s ) 

1 ^ nnx TTOX' e - nn (t-t')/L + e -nn(P-t+t')/L 

— > sin —— sin —-—-r- 

7 ra L L nil — q n ) 

71—1 V 7 


(4.20) 


The width is ultraviolet divergent and is regularized with the point-splitting method, 
where the primed coordinates are assumed to be e away from the non-primed coordinates, 
in the limit e —> 0 giving unmans] 


w 2 < k L/2) = —log^ + — log7](2r), (4.21) 

ZlI(J J-jq 7Y(J 

with Z/ 0 proportional to e. In the computation one finds essential the property 

OO n OO 

= “S log ^ _<?n ) = ~ lo s^( r ) > ( 4 - 22 ) 

71—1 71—1 
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where the function 0(r) is the Euler function, related to the Dedekind ^-function by 
4>{t) = 77 (r)g -1 / 24 . 


Equation (14.211) has two terms. When (3 L the leading term is the logarithm, and 
the logarithmic broadening is obtained. When L (3 the width becomes instead 

wHL/2) = ±. tog A. + A. L , ( 4, 23 ) 

and the leading term becomes linear. The fluctuations in the former limit are more inter¬ 
esting since are closer to the Wilson loop setup in terms of the limits used in holography. 


4.2.2 Curved backgrounds 


We still assume periodic time and a Euclidean signature on the worldsheet and introduce 
the notation 

Sx 2 (cr,t) = x(cqf). (4.24) 

In order to use the AdS metric, we are working in the low temperature limit LT -C 1. 
The equation per independent fluctuation reduces to 

d a ( ^7=t<9<xX(t t) \ - —— <9 t 2 x(a, t)= 0 , (4.25) 

\ \ r / £?00 

where gn and F are functions of a. We express everything in parameters x = a j Uo with 
x G [1, 00 ) and t = ta ^' 2 and rewrite the above differential equation as 

(PX')' ~rx = 0 , (4.26) 


where 


P 


9 00 9 11 


r = 


9nVF 


(4.27) 


V~F ’ Poo 

are functions of x only, the prime denotes differentiation with respect to x, the dot 
differentiation with respect to t and s is the total scaling dimension of the equation. The 
periodicity of the time direction (3 then becomes (3 = f3a~Q S ^ 2 . 

We are therefore looking for the Green function G that satisfies the equation 


( pG 1 ) 1 — rG = 5(i — t') S(x — x') 


(4.28) 


and obeying the jump condition 

- r (d(j3) - G(0)) = 5{x - x '), 


(4.29) 
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obtained by integrating eq. (14.281) over time. Next we define the expansions 


where 


XOM) = J ^c n (t)w n (x ) , 

n 

/ oo 

p(x' ,t)w n (x')r(x')dx' , 

n 

1 

^n(^) G n 7l n (x^j , iV n 


/i°° u n (x) 2 r(x)dx 
is a set of complete, orthonormal functions on [0, oo) satisfying 

/ OO 

dx w n (x)w m (x) r{x) = 8 nm , W re (x) w re (3/)r(a/) = 5(x -a/) . 


(4.30) 

(4.31) 


(4.32) 


(4.33) 


We have introduced the auxiliary source function r(x)p(x,t ) in order to compute the 
Green function. Inserting on the right hand side of (14.201) the source we obtain (summa¬ 
tion over n is implied) 


c n (pw' n )' - rc n w n = r(x)p n (t)w n (x) . 


(4.34) 


Requiring w n to solve 

( PK)' -fl 2 n rw n = 0 , 

eq. (14.341) is satisfied if 

C n G n C n Pn • 


(4.35) 


(4.36) 


The solution to the above, with the periodic boundary conditions in the time direction, 
is 


p Qn(t t') _|_ g (/3 (t t')) 


Cn(t)= (i- t') (-pn(O) dt' = / -——-- Pn(t')dt\ (4.37) 

Jo Jo 

where the Green function (t — t') depends only on the dimensionless time. Then, 
using eq. (14.301) we have 


x(M) = 

with full Green function 


/3 POO 



G(x, t; x', t) p(x', t')r(x') dx'dt' , 


(4.38) 


o J 1 


e -Q n (t-t') _|_ e -n n (j3-(t-t')) 

G(x, t; x\ t') = 22 -7 GTT-, - 7TT\ - w *( x ) w n(x ) ■ (4.39) 

„ 2G n (l -e-^P) 
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In terms of dimensionful time this is 


e -Ci„{t-t')a 0 s/2 _|_ g-Q n (P-{t-t'))o\ 


G(x, t; x ', t 7 ) = 


-s/2 




20 , 1 - e~ Qn ^o 


-s/2 


w n (x)w n (x') . 


(4.40) 


Notice that for the AdS metric analysis presented in the following section, s = —2, a 0 ~ 
1 /L and for the toy example presented in Appendix El s = +2, a 0 ~ L. the dimensionless 
time parameter is r = t/L and the Green function for both cases becomes 

G(x,t-,x',t , ) = Y^ - 2^(1 - e ~^/ L ) -(4-41) 

with w n satisfying eq. (14.351) . 


4.3 The width of the fluctuation 

Working with the metric elements gn = a 2 , g uu = 1/cr 2 we have 

x 2 

s = — 2 , r = ^ i , p = x 2 \fx 4 — 1 . (4.42) 

Writing (14.351) in the form of 

y" +-y'- G? n -y = o (4.43) 

P P 

the equation for the AdS Wilson Loop becomes 

y " + tV- 2 e) ) y ' + n " (x4 - 1)j/ = 0 ' <444> 

We need to find the Green function of this equation and perform the point splitting 
around x — 1. An analytic computation does not seem to be possible, but a numerical 
version with certain approximations may be done. 

To begin, we can solve numerically eq. (14.441) in order to specify the solutions w n up 
to normalization constants. Then using the formalism of the previous section yielding 
the Green function (14.391) . we can specify numerically the Green f un ction of (14.441) . The 
point-splitting is also done numerically and taken close to the value x — 1, where the 
result is integrated over all positive values of uj. Before listing the results we mention 
a few shortcomings of the numerics. Although we can go infinitely close to x = 1 we 
can not take the exact value of unity and therefore we are forced to introduce a cut-off 
x m \ n = 1 + e 2 . The integration in all positive values of u> should also start with a non-zero 
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value (u m i n = 6| < 1 and arbitrarily close to zero. Therefore, in the numerics we have 
introduced these two cut-offs even though the normalized result seems to depend weekly 
on their exact values as long as they are small enough. It should be finally mentioned 
that we work with large numbers, as expected in these limits, that could be the source 
of numerical uncertainties. 

Our numerics indicate that the leading contributions come from the time dependent 
part of the Green function. Moreover, we find that the width is linear in L/((3e i) to 
a good approximation. This is a term that exists also in the width of the toy model 
configuration presented in Appendix A and appears in equation (IA.17lh The cut-off 
independent term that we find to a good accuracy is the logarithm In (L//3), also found 
in the configuration of Appendix in flA.181) . Therefore, the total fitted expression reads 

” 2=ici ^“' :2in ? +F ( ei ' £2 ' £)■ <445) 

where F-> (ei, e 2 , j-) is an unknown function where its argument dimensions are normalized 
by the radius R of the AdS space where needed, and ci, c 2 > 0 are constants that are 
specified numerically. By adding higher powers of L//3 to the fit it can be seen that the 
values of the coefficients ci, c 2 are stable. We remark that our results are valid in the low 
temperature limit, which is translated to L <C /3. Moreover the fitted expression should 
be taken with caution due to the peculiarities of the numerics we have mentioned and 
the presence of the undetermined function F. We hope to return to this subject in the 
future with more advanced results. 


5 Conclusions 

In this paper we have studied the properties of the flux tube connecting static quarks 
and anti-quarks in several cases. We have used two alternative approaches. 

Firstly, working with generic confining theories we have calculated the flux generated 
by a heavy quark pair through a probe Wilson loop. This is achieved by placing the 
measuring small Wilson loop in the middle of the heavy bound state and at a distance 
above it. The resulting computation in the large N strong coupling limit turns out to 
be a connected minimal surface with boundaries the two Wilson loops. The width of the 
tube in the middle of the pair can be found by specifying its minimal area. Using the 
properties of the string solutions and the minimal surfaces in confining backgrounds, we 
have argued that the logarithmic broadening of the flux tube is a property of all confining 


19 






theories in the context of the gauge/gravity duality. Our analysis is supported by numer¬ 
ical computations, where the system of equations is solved without any approximations. 
Notice that our analysis is not applicable to conformal backgrounds and it will be sur¬ 
prising to us if the same properties will be observed in backgrounds where confinement 
is not present. The reason is that the minimal surfaces will have different shape when 
confinement is not present. 

In the large quark separation limit the rectangular Wilson loops can be replaced by 
circular ones and this is not expected to affect the type of broadening. If one would like 
to repeat our analysis in non-confining theories with the presence of black hole at finite 
temperature, this assumption needs to checked. It would be also challenging to build 
consistently the tubes in the Lifshitz-like theories and in generic anisotropic theories at 
finite temperature, like for example in Pi, where the calculation of several observables 
involving the interaction between the heavy quark pair has been made [HD E21E3. Those 
backgrounds may be also considered as modified hard wall models to secure co nfi nement. 
The anisotropy of space will impose a non-trivial shape on the minimal connecting surface 
and it would be interesting to see to what extend the logarithmic broadening of the 
width will hold. Notice that in the anisotropic theories have been found several universal 
properties of the isotropic theories to be violated (HUES]. 

As an alternative approach to understand the dynamics of confining strings we have 
considered fluctuations of strings in particular limits and in curved backgrounds. The 
strings have fixed end-points, and their time evolution leads to fluctuations in space 
time. We started with the time-independent fluctuations of strings with fixed end-points 
on the boundary. Then time-dependent fluctuations have been considered in curved 
backgrounds. The Green function of a differential equation with particular boundary 
conditions needs to be found in order to specify the width of the fluctuations. The 
setback is that the differential equations of the fluctuations can not be solved analytically. 
We have employed a point-splitting regularization of the Green function and carried 
out numerically the computation in the simplest case of the curved AdS background. 
Although the dual theory is non-confining we computed an analogous quantity to the 
width of the heavy meson. There are difficulties in finding the two-point function even 
numerically, mostly due to the singularities which generate certain divergencies that need 
to be isolated. In the numerical fitting we were able to see that for small distances the 
quantity we have defined as width seems to have at least a logarithmically reducing term 
with the distance of the string, and that would be interesting to see if this is a generic 
property of non-confining theories. Due to some peculiarities of the numerics we plan 
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to return to this result in the future for a more accurate treatment. As a next step one 
would like to do the full analysis of fluctuations in order to compute the width of flux 
tubes in confining theories. 

An extension of our work would be to study the broadening of the tube between 
more involved heavy bound states formed in gauge/gravity duality than the quarkonia. 
The strongly coupled bound state formed by k interacting heavy QQ, called k-strings is 
expected to have the logarithmic broadening for any confining theory. More demanding 
may be the holographic computation of the baryon effective width, which in effective 
string models using the Y ansatz has been found to grow logarithmically [13j. Holo¬ 
graphically the calculation maybe done with the baryon vertex string [22 E3, which 
already has ben found to have certain universal features among the different theories 
[38j. Moreover, it would be very interesting to relate the entropy grow associated with 
the heavy quark-antiquark bound state with the properties of the flux tube along the 
lines of m m El- 
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A An analytically solvable toy model 

In this section we use the mathematical formalism developed in the main text, to a toy 
model using a space-time string that leads to analytic results. The motivation is just 
to present an example where the Green function problem can be solved and its point 
splitting regularization can be carried out till the end analytically. So we present a toy 
model of a string embedded in the curved space, where the choice of its metric has been 
motivated by the fluctuation equations (14.251) . in order to obtain analytical solutions. 
The metric of the space is ds 2 = \/l + u(dx 2 + dt 2 ) + (vTTa) _1 tf. The analysis of the 
string is the same as in section 14.11 with the difference that the u direction has different 
scalings in the curved space where the string lives. In order to keep the length of the 
string finite in this case we need to introduce a large upper bound M on a at the direction 
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u, where we apply the boundary conditions. In this simple space we find the fluctuations 
of the string. We point out that this is only a toy working model, which however allows to 
develop the analytic methodology, and exhibits some interesting features of the defined 
’’width”, which we compute. The turning point of the string is related to its length via 
eq. (14.91b in this case by 

L 


a ° 2W_i(—L/4M) ’ 

where W- i(x) is the generalized Lambert function. The equation of the transverse fluc¬ 
tuations that need to be solved is (14.351) in the form 


y"+ P -y'-tf n -y = 0, 


with 


s = 2 , r = 


P 


x 


P 


P 


= x 2 \/x 2 — 1 . 


(A.2) 


(A.3) 


\/x 2 — 1 

We observe that the singularity at a = cr 0 , that is at y = 1. The solution to eq. (IA.2I) 
can be found if we set y = (1 — a; 2 ) i/4 /i(.i:)/.i:. Then we obtain for the unknown function 
h 


1 


(1 — or) h" — 2 xh' + 


' 2 3 1 

(V -I- 

n 4 4(1 — x 2 )_ 


h\ =0. 


(A.4) 


x(l — x 2 ) 3 / 4 

which is the associated Legendre equation with solutions the associated Legendre func¬ 
tions P^j 2 and QxJ- where A n = —1/2 +K n and K n = -y/l + fl 2 . Therefore the differential 
equation eq. (IA.2D has the two independent sets of solutions 


— (1 — x2 ) l ^ A d i /2 _ — x2 ) 1 ^rP/ 2 

U l,n ( x ) — ^ P\ n > u 2,n( x ) ~ ^ Q\ n > 


(A.5) 


which can be written in the simple algebraic form 


P A 1/2 (x) = _ (x 2 — 1) l P l x T \/x 2 — 1 

An V^TT 


An + 1/2 


+ 


[x + — lj 


An “I - 1/2 


Q\?{ x ) = 


Ut 2 _ 11-1/4. 


1 


[x + \J~xP — lj 


An + 1/2 


(A.6) 


By making the change of variables x = cosh z , with z G [0, oo) our variable x is indeed 
G [1, oo), we get the compact expressions 


Ul, n (z) = 


2 cosh(W n 2 ;) 
7T cosh ^ 


U2,n( z ) = 


7r e 


—K n z 


2 cosh 2 ; 


(A.7) 
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The first set of solutions Ui jH (z) is ruled out using the boundary conditions and require¬ 
ment of normalizability and we end up to a Green function of the form 


G(z,t',z',t') = 


e -Cl n (t-t')/a 0 _|_ e -Cl n (/3-(t-t'))/ao e ~K n z e ~K n z' 

^ 2f2 n (l — e~ flnl3 / rTo ) cosh 2 cosh z 1 


(A.8) 


where we have used the generic function (14.3911 and /V 2 = K n by the normalization 
condition (14.3211 . 

We continue with the point splitting regularization of G by setting t 1 = t + e, z' = 
z + e /00 and adding and subtracting in the numerator a term in order to split G into a 
sum of a /3-dependent and a /3-independent part. This leads to 


= W lo + w l 


(A.9) 


with 




= £ 


2ft, 


—e v CT o J -\- • • • 


2 i V n 


K n e Kn °o cosh(f2 n ^) 
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1 - e n CT o 


+ • • • 5 


(A.10) 


where we have kept the leading terms from the expansion of the cosh e dependent terms. 
In the limit of a continuous spectrum we take fl n —> u> and then to leading order we 
have 
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00 , VT 

auo- 


2uj 


—-e J 
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/ d u -;- r~ ■ ( A - n ) 
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1 — e CT o 


We treat separately the two integrals, focusing first on w^. With the change of variables 
h = \/l + a; 2 — co we obtain 


2 f 1 J u (h 2 + 1) - e -h 

= / ah -—e CT o 


(A. 12) 


/ 0 4h 2 (/i 2 — 1) 

where integrating the above and expanding in e, we get elementary functions and the 
exponential integral function. In the small e limit we arrive at 


= - (l — log 2 — lim log h] —— ( - + 4log2 + 6 lim log h ) . 

2 V h —^0 / (Tq V 8 h —^0 


(A.13) 


Next we turn to w where we notice that the dominant contributions come from the 
uo 0 regime. By introducing a cut-off r 0 in the integral, the expression in flA.111) 
can be written in the form 


w l ~ 


duo 


\/l + —— x/l+cd 2 \ A — —(n+l)uo 

— ^ ^ rr n \ / 


-e CT o 


UJ 


£ 

n=0 


e a o 


(A.14) 
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Figure 6: The quantity w 2 in terms of L, where for small distances where our approximation 
is valid is reducing. For the plot we have fixed the normalized quantities M — 600, /? = 400, 
which gives ctq «M. 


where the dominant contributions come from the uj 0 regime. We expand in e and we 
integrate to obtain 


(X) 

w l = Y. 

71=1 


r o, 


f3nr 


'o 


do J 2f3 2 n 2 
that can be written in the more convenient form 


+ ^C(2) 



2/3 2 


1 - — 
Vq 


1 - 


v 0 


(A.15) 


(A.16) 


where £(2) is the zeta function. Using the Euler-Maclaurin formula, we get for the leading 
contributions 


w p = 


+ 7 - 1 + In — + lnr + ^ 
pr (j 0 2 p 2 


C(2) 


1 - 


vo 


(A.17) 


with 7 the Euler constant. Isolating the divergences coming from the cut-offs we find 
from (IA.13D (IA.17I) the total expression for the ’width’ 


2 1 Vq 

w = Ci — in — + 


7T 2 ok 


f3 12 f3 2 (j 0 


+ — ( c 2 + In ^ - 


7 T 2 ok 


13 12 (3 2 


(A.18) 


with Ci and c 2 constants that can be easily specified and o 0 given by flA.ll) . In Figure El 
we observe that at the small L and cy < M the quantity we treat as width is reduced in 
this case with the increase of length. 
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